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ABSTRACT 

The recently discovered population of ultra-faint extended line emitters, with fluxes 
of a few times 10~^^ erg s~^ cm~^ at z ^ 3, can account for the majority of the 
incidence rate of Damped Lya systems (DLAs) at this redshift if the line emis- 
sion is interpreted as Lyg. We show here that a model similar to that proposed by 
iHaehnelt. Steinmetz. fc RauchI (120001 ) , which explains the incidence rate and kinemat- 
ics of DLAs in the context of ACDM models for structure formation, also reproduces 
the size distribution of the new population of faint Lya emitters for plausible pa- 
rameters. This lends further support to the interpretation of the emission as Lya, as 
well as the identification of the emitters with the hitherto elusive population of DLA 
host galaxies. The observed incidence rate of DLAs together with the observed space 
density and size distribution of the emitters suggest a duty cycle of ^ 0.2 — 0.4 for 
the Lya emission from DLA host galaxies. We further show that Lya cooling is ex- 
pected to contribute little to the Lya emission for the majority of emitters. This leaves 
centrally concentrated star formation at a rate of a few tenths M0yr~^, surrounded 
by extended Lya halos with radii up to 30-50 kpc, as the most plausible explanation 
for the origin of the emission. Both the luminosity function of Lya emission and the 
velocity width distribution of low ionization absorption require that galaxies inside 
Dark Matter (DM) halos with virial velocities < 50 — 70 km s~^ contribute little to 
the incidence rate of DLAs at z ~ 3, suggesting that energy and momentum input due 
to star formation efficiently removes gas from these halos. Galaxies with DM halos 
with virial velocities of 100 — 150 km s~^ appear to account for the majority of DLA 
host galaxies. DLA host galaxies at z 3 should thus become the building blocks of 
typical present-day galaxies like our Milky Way. 
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1 INTRODUCTION 

Quasar absorption spectra provide excellent probes of 
the distribution of baryons in the high-redshift Universe. 
Damped Lyman Alpha systems (DLAs, historically de- 
fined as having a neutral hydrogen column density A'hi > 

2 X 10^°cm~'^) are particularly useful as they are likely to 
play an important role as a reservoir of gas for the for- 
mation of stars and galaxies at high redshift. They domi- 
nate the neutral gas content of the Universe between z ~ 
— 5, and at z ~ 3.0 — 4.5 their neutral gas content is 
cornparable to visible stellar mass in present-day galax- 
ies (jW olfc 1986; Storric-Lombar di. McMahon. fc Irwinll996l : 
IStorrie- Lombard! fc Wolfell2000l ) . DLAs thus form an impor- 
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tant link between primordial plasma and the stellar struc- 
tures that form from it. 

In spite of observati o ns of over 1000 DLAs, 
IWolfe. Gawiser. fc Prochaskal (|2005D still conclude that 
the question "what is a damped Lya system?" has not yet 
been answered conclusively. One of the reasons for this is 
that absorption spectra provide only indirect information 
in velocity space, and only probe the gas along one line- 
of-sight through the galaxy. The Lya absorption feature 
itself provides no information about the velocity structure 
of the DLA, because of the large optical depth even in 
the damping wings. This has led observers to look at the 
associated low ionization metal absorption features. Low 
ionization species like Sill, CII and Fe l are believed to be 
good tracers of the neutral gas in DLAs. iProchaska fc Wolfel 
(| 19971 ) developed the velocity width distribution of metal 
absorption features into an important diagnostic tool for 
DLAs. The much lower absorption optical depth allows us 
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to extract detailed kinematical information for the gas in 
DLAs. The absorption profiles are normally clumpy and 
asymmetric, with the strongest absorption feature often 
occurring at one edge of the profile. Velocity widths range 



from 30 km s up to several hundred km s 



Note, 



however, that there is a distinctive lack of very narrow 
absorption profiles. 

Kinematical models aiming to reprod uce the veloc- 
ity width data fall into two categories. IWolfd (|l986l ) 
suggested a close connection between DLAs and disks 
of present-day spiral galaxies. iProchaska fc Wolfd (Il997l ) 
modelled DLAs as thick rotating disks with a rotation 
speed 200 km s~^) typical of present-day galaxies. 
iHaehnelt. Steinmetz. fc RauchI (|l998f ) challenged this inter- 
pretation and demonstrated that it is not unique. The 
merging of proto-galactic clumps expected in Cold Dark 
Matter (CDM)-like models for structure formation can ex- 
plain the shape of the profiles equally well. Galacti c winds 
are also likely to play an important role (e.g. iFabianI 
1 19991 ). They further showed that for the same virial ve- 
locity width, merging clumps produce significantly larger 
velocity width and argued that the latter interpretation 
is favoured by the observed velocity width distribution in 
the context of the CDM paradigm Ib r structure formation 
jHaehneh. Steinmetz. fc Rauchll200ol ) . 

The level of enrichment with metals provides an- 
other ini£ortant_clue as^^ — see 
IPettinil (|2004l . l2006l l: IWolfe. Gawiser. fc Prochaskal (|2005l l 
for reviews. Significant metal absorption is found in all 
DLAs, though as a population the y are metal poor . 
DLAs are also relatively du st free (jEllison et al.l l200ll : 
I Wild. Hewett. fc Pettinill2006l ). Initially, the low metalhcity 
was used to argue that DLA host galaxies are the chemically 
unevolved but otherwise very similar counter-parts of typi- 
cal present -day spiral galaxies llWo lfc 1986). In the model of 
iHaehneltT Steinmetz. fc RauchI ( 1998|) , DLAs instead prefer- 
entially probe the outer parts of much less massive galaxies, 
many of which end up as building blocks of typical present- 
day galaxies that form by hierarchical mer ging in CDM- 
like m odels for structure formation. Recentlv. lPontzen et al.l 
l|2008h demonstrated that such a model fits the observed 
metalhcity of DLAs at z ~ 3 very well. 

Models of DLAs invariably predict more than quasar 
absorption features can in principle test. As already men- 
tioned, absorption is measured along a single line of sight 
and thus is a one-dimensional probe of the properties of 
the DLA. To explore the spatial extent and structure of 
the DLA, we need to observe DLAs in emission. Attempts 
to do this have focussed on both line and continuum emis- 
sion. Lya emission holds great potential in this respect. Star- 
formation, cooling radiation and fluorescent re-emission of 
the meta-galactic UV background are all expected to con- 
tribute at different levels. In addition, stellar continuum 
emission from the newly formed stars should also be bright 
enough to be detectable. 

Many observers have attempted to find the galaxy 
counterparts of DLAs at high redshift in emission by 
searching adj acent to quasar sightlines with known absorp- 
tion sy stems llFvnbo. M0ller fc Warrenll 19991: iBunker et alj 



19991: iKulkarni et al.l |2000|, l200ll: IWarren et al.l l200ll: 



tracted to stud y the light of th e galax y, which is very faint 
in comparison. iKulkarni et al.l l|2000l ) and others caution 
of the possibility that a given emission feature is a Point 
Source Function (PSF) artefact rather than a real source. 
IChristensen et ahl l|2007h report that, for 2 > 2, six DLA 
galaxies have been confirmed through spectroscopic obser- 
vation of Lyof emission, w ith other techniqu es producing a 
few additional candidates. IChristensen et ahl added another 
six LyQ emission candidates to this group. A quantitative 
statistical interpretation of the many (largely unsuccessful) 
searches is difficult if not impossible, but the rather low suc- 
cess rate appears to be consistent with their interpretation 
as galaxies of rather low mass and star formation rate. 



iRauch et al.l l|2008l ) reported the results of a long-slit 
search for low surface brightness Lya emitters at redshift 
z ~ 3, which reached fiux levels that are about a 10 times 
lower than previous Lya surveys at this redshift. They found 
27 faint line emitters, many of which are extended in wave- 
length and real space. They argue that the majority of the 
emitters are likely to be Lya (rather than \ow-z interlop- 
ers) , powered by a central region of star formation and pro- 
cessed by radiative transfer through surrounding neutral hy- 
drogen. (Pauch et al.l l|2008l ) note that the incidence rate in- 
ferred from the space density and the size distribution of the 
emitters is similar to that of DLAs, and suggest that they 
are the host population of DLAs and high column density 
Lyman Limit Syste ms (LLS). If this is i ndeed the case then 
the observations of iRauch et al.1 (|2008h give us — for the 
first time — the size distrib ution and s pace density of DLA 
host galaxies (Figure 19 of iRauch eni.l.ii2008i '). 



Models predicting the velocity width and size 
distribution of DLAs have been constructed based 
on the observed luminosity function of galaxies 
iFvnbo. M0ller. fc Warrenlll999l : iFvnbo et al1l2008l) and on 
the Press-Schechter formalism ( Press fc S chechtcr 197'|) 
in con junc tion with numerical simulation s (Gardner ct al.l 
| l997allbl: IHaehnelt. Steinmetz. fc RauchI Il998l. 



Christensen et al. I I2OO7I ). This is a difficult task, as the light 



of the extremely bright quasar must be accurately sub- 



Nag amine. Springel. fc HernauistI l2004l : iNagamine et al.l 
2007). More recently numerical simulations have attempted 
t o model the entire DLA population self-co nsistently 
(|Razoumov et al.ll2006l . l2007l : IPontzen et al]l2008l l. 



In this paper, we revisit models for the absorption prop- 
erties of DLAs in t he light of the size distribution data of 
iRauch et al.l l|2008l ). improved data on the velocity width 
distribution from metal absorption lines (which has pre- 
sented a challen g e to purely numerical simulations — see 
Razoumov et al.| (|2007l )) and numerical simulations with in- 
creased resolution, box size and sophistication incorporating 
the additional physics of radiative transfer, gas chemistry 
and star formation. 



Throughout this work we use cosmolog- 
ical paramete rs of the 5 year WMAP data 
(jHinshaw et al.l l2008h : {h,0,M,Qb,^A,o-s,n) = 
(0.701, 0.279, 0.046, 0.721, 0.817, 0.96). 
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2 A JOINT MODEL FOR THE KINEMATICAL 
PROPERTIES OF DLAS AND THE 
CUMULATIVE SIZE DISTRIBUTION OF 
THE FAINT Lya EMITTERS 

2.1 The Haehnelt et al. model 



iRauch et alJ (|2008h argued that their population of faint 
LyQ emitters is the same or has at least a large over- 
lap with that of DLA/LLs host galaxies. They further 
pointed out that the space density and sizes should 
agree w ell with those pred icted by the DLA model of 
iHachncl t. Steinmetz. fc Ranc h (1998), which models DLAs 
in the context of CDM models of structure formation. 

We here revisit and update this model to investigate 
whether it can explain the properties of DLAs and the new 
population of faint emitters, assuming that these are the 
same objects. We start with summarizing the salient prop- 
erties of the model. As discussed in the introduction, ab ini- 
tio numerical simulation of the gas at the centre of galaxies, 
where complex non-linear gas physics including star forma- 
tion and the associated feedback are important, is still very 
challenging. The Haehnelt et al. model therefore takes a hy- 
brid approach, using a combination of Press-Schechter for- 
malism and results from numerical simulations to model the 
kinematic properties of DLAs (see [Johansson fc Efstathioul 
l|2006ll for a semi-analytical model of DLAs that explicitly 
models feedback). 

The model uses the space density of DM halos as a 
proxy for the space density of DLA host galaxies. With the 
refinement to the Press -Schechter formalism introduced by 
ISheth fc TormenI (|2002l ). the number of dark matter halos 
per unit comoving volume at redshift z with mass (baryonic 
-I- CDM) in the interval (M, M + AM) can be estimated as. 



2 po Av' 



nM (M, .) dM = A ^1 + — J ^ - ^ — exp ^- - J AM, 

(1) 

where um is the rms fluctuation amplitude of the cos- 
mic density field in spheres containing mass M, po is the 
present cosmic matter (baryonic -I- CDM) density, i/' = ^/ai^, 
V = 5c/\D(z)aM]- D(z) is the growth factor at redshift 
z (|Carroll. Press, fc Turneilll992l ). Sc = 1.686, a = 0.707, 
A ~ 0.322 and q = 0.3. We have used the fitting formula in 
lEisenstein fc Hul |l99i) to calculate the matter power spec- 
trum. 

We will be interested in the kinematic properties of 
DLAs, for which the virial velocity is a more convenient 
quantity to characterize the DM h alo than the mass . The 
two are related as follows (e.g. Mail er fc Bullockl 120041 '), 
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where Ay is the overdensity of the halo (see 
iBrvan fc Normanl |l99i). 

We use a simple power-law relation between virial ve- 
locity of DM halo and absorption cross section 



200 km s- 



(3) 



where (3 and ro are parameters. The value of f3 has been 
the source of some controversy, due to numerical simu- 
lations still finding it challenging to reliably model the 
spatial distribution of th e gas in the h i gh-den sity region 
probed by DLAs and LLs. [Gardner et al.l l|l997al ') originally 
favoured a value of /? = 2.94 (at z = 3, for a ACDM cos- 
mology), but their later work revised this to (3 = 1.569. 
[Prochaska & Wolfe (2001) pointed out that this low value 
of is incompatible with the ob s erved DLA velocity widths. 
iHaehnelt. Steinmetz. fc RauchI (|200G| ) found that the ob- 
served velocity width distribution of metal-lines is repro- 
duced well with a value of /3 ~ 2.5, which we also use here. 

Haehnelt et al. further found that they needed to in- 
troduce a lower cut-off in virial velocity in order to fit 
the velocity width distribution of low-ionization absorption 
systems. Otherwise their model predicted too many very 
narrow metal absorption systems, which are not observed. 
They therefore assumed that DM halos with virial velocities 
smaller than a minimuin velocity fmin do not host DLAs. 

Here we will slightly relax this assumption and model 
the suppression of the cross-section of DLAs in halos with 
small circular velocities as a gentler exponential decline. 
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We consider a range of values for the parameters Vcfi and a. 
Note that a sharp cut-off corresponds to a = oo. 

It remains to choose the parame ter ro. We will follow 
IHaehnelt. Steinmetz. fc RauchI l|2000l ) by fixing ro so that 
the overall rate of incidence of absorbers per unit redshift 
dTV /Az agrees with the observational value, 
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We take the value of dTV /Az = 0.24 at z = 
3 , which is consistent with AJ\f /A X = 0.067 
(iProchaska. Herbert-Fort, fc Wolfd l2005l : iPeroux et all 
l2005h . The ratio of proper distance interval to redshift 
interval is given by. 
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and the so-called absorption distance, X, is defined by. 
Ho 
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(6) 



(7) 



The resulting values of ro are given in the caption of Figure 

m 

Further, it is useful to define the line density distri- 
bution (N), such that the number of intersections (A^JV) 
of a random line of sight with systems that have mass in 
the interval {M,M + AM), located in the redshift interval 

(z, z + Az) is. 



dW = N(M, z) AM Az 
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Az 



The velocity width distribution of 
low-ionization absorption 



(8) 
(9) 



To calculate the velocity width distribution, Haehnelt et al. 
used a conditional probability distribution: p{v^\vc) duw is 
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Figure 1. Top Panels: The velocity width distri bution l(v„,X) of the associa t ed low -ionization metal absorption of DLAs. The black 
crosses show the observational data compiled in IWolfe. Gawiser. fc Prochaskal l|2005l V Bottom Panels: The cumulative velocity width 
distrib ution v^,X) of the a ssociated low-ionization metal absorption of DLAs. The thick solid black cu rve shows the observationa l 
data of lWolfe. Gawiser. &: Proc haska (2005). The thin, solid black curve shows the older observational data of lProchaska &: Wolfd lll997t ). 
Left Panels: The velocity width distribution of our model for an exponential suppression at small virial velocities of the cross section for 
damped Lya absorption as given by Equation (jjjl with a = 3 and Vcfi = 30, 50, 70 km s~^, respectively. The cross section normalization 
for these models as given by Equation JSj is characterised by ro = 18.5, 23.5, 28.4 kpc. Right Panels: The velocity width distribution 
for Vcfi = 50 km and a range of values of a = 1, 2, 3, oo. For these models, ro = 21.9, 23.6, 23.9, 23.5 kpc. 



the probability that a DLA in a halo with circular velocity 
Vc has a velocity width in the interval (fwi^w -I- dww). On 
the basis of their numerical simulations, they assumed it to 
be a function of the ratio Uw/wc- The number of systems 
along a random line of sight in the interval {X, X + AX) 
with velocity width greater than Ww is given by 



l{>v^,X) = 



p{v^\vc{M)) 



A?M 
AX AM 



AM 



Av^ 



(10) 

We will u se here the conditional probab i hty d istribution 
as given in lHaehnelt. Steinmetz. fc RauchI l|200d . Figure 1). 
The distribution peaks at v^/Vc ~ 0.6, dropping to zero 
below Vv,/vc ~ 0.1 and above v„/vc ~ 2. The numerical 
simulations on which the distribution is based did not con- 
tain star formation feedback. Simulations have become more 
sophisticated since then. We have therefore compared the 
p{vvi\vc) distribution used here with that from the simula- 



tions of iPont zen et akl l|2008l . and private communication), 
which incorporate the effects of star formation and super- 
novae on the kinematics and spatial distribution of the gas in 
a simple manner. The differences in p{v^\vc) axe small. The 
resulting difference in Z(> Vv,,X) is also small, at most 10%. 
Unless there is a fortuitous cancellation of different effects, 
this suggests that star formation in the simulations has a 
small effect on p{v^\vc). This is somewhat surprising, given 
the significant differences in resolution, cosmological volume 
and additional physics, albeit reassuring for our modelling. 
Note, however, that the simulations still fail to produce re- 
alistic galactic winds, probably due to the rather simplistic 
fashion in which stellar feedback is incorporated. 

Figure [T] compares the predicted velocity 

width distribution with the obs ervational data of 

Wolfe. Gawiser. fc Prochaskal (|2005D . The data of 

Prochaska fc Wolfd is also shown in the 

bottom panels and is reasonably well fit with 
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Figure 2. The cumulative incidenc e rate inferred from the size distribution of the population of faint Lyo emitters. The solid black 
curves are the observational data of lRauch et al Left: Cumulative incidence rate in our model for o = 3 and a range of values 

of ''c.o = 30, 50, 70 km s~^. Right: The cumulative incidence rate in our model for v^fl = 50 km and a range of values of the duty 
cycle /d = 1, 0.5, 0.2. 



fc.o = 50 km s ^ . This is very s i milar to what was found by 
iHaehnelt. Steinmetz. fc RauchI (l2000ll. The new compila- 
tion o f velocity width data by Wolfe. Gawiser. fc Prochaskal 
(|2005l l extends to significantly larger velocities and appears 
to require a somewhat larger value of Wc,o = 70 km s~ . 
The apparent lack of neutral gas in small dark matter halos 
can be plausibly attributed to the feedback effects of star 
formation and/or photo-heating due to the meta-galactic 
UV background, fn most numerical simulations and models 
of galaxy formation, these feedback mainly affects halos 
with somewhat smaller virial velocities than this. So either 
halos with small virial velocities have larger velocity widths 
than we have assumed here (i.e. p(«w|i'c) is different), 
or else stellar feedback in halos with virial velocities of 
Vc,o ~ 50 km s~^ is more efficient than generally assumed. 
The cumulative velocity width distribution is not very 
sensitive to the shape of the cut-off for q ^ 2. 

The hi gh velocity tail of the velocity dist ribution as 
compiled bv lWolfe. Gawiser. fc Prochaskal (120051 ) has proven 
difficult to reproduce with numerical simulations, which at- 
tempt to model the spatial distribution of the neutral gas in 
DLAs self-consistently rather than assume a scaling of the 
absorption cross section with the virial velocity of DM halos. 
Generally, the simulations fail to produce a sufficient num- 
ber of absorption systems with velocity widths as wide as 
observed. This is normally attributed by the authors of these 
studies to the fact that momentum and energy input into the 
gas due to star formation may not have been modelled with 
suffic i ent sophistication (iNagamine. Springel, fc Hernguistl 
20041: iNagamine et al.l l2007l : iPazoumov et all |2006| , I2OO7F 



Pontzen et al. 



2008h . 



2.3 The cumulative size distribution 

We now ask whether our model, which successfully fits the 
velocity width distribution of the associated low-ionization 
metal absorption of DLAs, can also reproduce the cumula- 
tive size distribution of the emission regions of the new pop- 
ulation of faint LyQ emitters, as shown by black solid curve 



in Figure (|Rauch et all |2008| . Figure 19, Ifl corrected). 
It asymptotes to dM /Az = 0.23, which is very similar to 
the observed incidence rate of DLAs at the same redshift. 
To calculate the cumulative incidence rate for the DLA host 
galaxies in our model, we must relate the size of the emis- 
sion region to the mass of the corresponding DM halo. We 
do this simply by assuming that the absorption cross-section 
is related to radius as if it was a sphere. The cumulative in- 
cidence rate is then given by 

AN , 



Az 



i> r,z) 



N{M',z) AM'. 



(11) 



M(r) 



The result is showrQ in the left panel of Figure (2] The pre- 
dicted sizes are consistently smaller than the observed sizes 
by a factor of between 1.5 to 3, suggesting that the emission 
regions of the Lya emitters are larger than the cross section 
for damped absorption. Even though the integrated inferred 
incidence rate is similar to that of DLAs, the emission re- 
gions of the Lyo! emitters can thus not be identical to the 
regions responsible for damped Lya absorption. Neverthe- 
less, the two can be closely related. 

To demonstrate this we now explore a simple model 
where only a fraction /d of the halos are emitting Lya radi- 
ation above the detection threshold at any one time. To keep 
the total inferred incidence rate fixed, we allow the cross sec- 
tion of the individual Lya emission regions to be larger than 
those for damped Lya absorption by a corresponding factor 
/J"^. The Lya emitters show typical signs of radiative trans- 
fer effec ts, like large veloc ity widths and asymmetric line 
profiles (|Rauch et al.ll2008l . Section 6.3), and should thus be 
optically thick to Lya radiation. However, we have no good 
handle here on the actual column density required to scat- 
ter emitted Lya at large radii effectively. It may well be 
lower than that required for damped absorption, in which 
case it is very plausible that the region for Lya emission ex- 
tends beyond that for damped Lya absorption. In the right 



^ We do not show the impact of altering a here. The effect is 
small, especially beyond a radius of lOkpc. 
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Figure 3. The cumu l ative luminosity function of the Lyo emitters. The black solid curve in both panels is the data taken from 
Fig. 9 of iRauch et al. The grey, shaded region is an estimate of the 1-cr confidence interval. Left: The blue solid curve and 

the dashed red curves are for optimistic models for the Lya emission due to Lya cooling radiation as described in the text. Right: A 
simple model for the Lyo emission due to star formation as described in the text. The relevant parameters for = (1, 0.5, 0.2) are 
(io'^minl km s-l]) = (4 X 10-*!, 75), (8 x 1041,60), (1.8 x 10-*2,45) respectively. 



panel of Figure [5] the cumulative incidence rate is shown 
for three values of the duty cycle fd ~ (1, 0.5, 0.2) with 
Wc,o = 50 km and a = 3. Our simple assumption of 
a duty cycle for the Lya emission reconciles the cumulative 
incidence rate predicted by our model with the observed dis- 
tribution for sizes below 35 kpc for fd = 0.2 {fd = 0.4 for 
Vc,o = 70 km s^i). The sudden drop in the observed distri- 
bution at r ~ 35kpc is likely attributable to the following 
two effects. The first is the surface brightness limit of the 
instrument. Light from sources with large radii may be too 
diffuse to be detected. The second is the eff'ect of searching 
a small survey volume. Large systems are rarer, so there is a 
limit to the size of sources that can be expected to be found 
within the rather small survey volume. 



3 MODELLING THE LUMINOSITY 

FUNCTION OF THE Lya EMITTERS 

3.1 The contribution of Lya cooling radiation 

I Rauch et al.l (|2008l ) considered a number of astrophysical 
origins for the Lya emission that they observe. They con- 
clude that the most likely mechanism for producing Lya 
photons is star formation. We here consider in more detail 
the argument that cooling radiation is unlikely to be the 
domina nt source of Lya for the observe d emitters. 

iDiikstra. Haiman. fc Spaanj (|2006l ) derive the following 
formula for the Lya luminosity {L^ya) of a collapsing pro- 
togalaxy due to cooling radiation, assuming that the gravi- 
tational binding energy is radiated as Lya on a dynamical 
time scale, 

5/3 



LLyc, = 5.8 X 10" 



Mu 



10" 

1.2 



1 -I- 2vir 



5/2 



2.5 



(12) 



erg s 



the halo, z^h is the redshift at which the system virialises, 
and the bulk velocity of the infalling material, Vhuik{r), is 
parameterised by Uamp and ad as a power law, fbuik(7') = 
Vs,inp{r /rvir)"'^ , where rvir is the virial radius0 We will set 
ZviT ~ 3, i^amp ~ Vc and consider the lower limit of th e 
range of ad discussed in lDiikstra. Haiman. fc Spaansll|2006h . 
namely ad = —0.5, so that we have an upper limit on I/Lya- 
Equation (|12|l gives a relation between the mass of a 
halo and the luminosity due to Lya cooling, which we can 
combine with the Press- Schechter formalism of Section [5] to 
predict the number of DLAs per unit comoving volume with 
Lya luminosity greater than some L^ya, 



-f'Lyai 2) — 



nM{M',z) dM' 



(13) 



v) 



There is, however, an inconsistency here. As already 
mentioned. Equation (|12|) assumes that the cooling radiation 
will be emitted over the dynamical time tdyn of the halo, 

^ - ( \ ^ 353 Myrs. (14) 

The proper time corresponding to 
[3.75, 2.67] of the Lya survey is = 



the redshift interval 
789 Myrs. Thus, all 
the halos cannot radiate at the luminosity L^yd given by 
Equation (|12|l for the entire redshift interval over which ob- 
servations were taken. There are two ways to make the model 
consistent. The first is to impose a duty cycle, so that at any 
one time, only a fraction fd ~ tdyn/iz ~ 0.45 of the halos 
will be radiating. The second is to reduce the luminosity of 
each emitte]0, so that the gravitational energy of the halo 
is radiated over 789 Myrs. 

The left panel of Figure [3] compares our prediction of 
n{> I/Lyai z) with the observed luminosity function as given 



where Mtot is the total (dark matter -I- baryons) mass of 



^ See Equation (10) of lDiikstra. Haiman. fc SpaansI 1I2OO6I) for a 
correction to this formula when Od < a^nd r is small. 
It could, of course, be a combination of these two effects. 



Emission and absorption properties of DLAs 7 



JO 



-1.5- 



(km s ) 
50 100 



500 




Razoumov et £ 
Pontzen et al. 
= 30km/s 



CO 



= 50km/s 
= 70km/s 



10 11 
log,„(M/M ) 

^10 ^ sun' 



13 



(km s ) 
50 100 




10 11 
'°9,„(M/M_) 



Figure 4. The contribution of different mass ranges and virial velocity ranges to the incidence rate of DLAs (d'^A/'/dX/d logj^Q M). 
The black squares in both panels show the results from the numerical model of iRazoumov et alj | |2007| ) f or DLAs found with in DM 
halos, excluding intergalactic DLAs. The grey symbols in both panels show the results of the simulations of I Pontzen et al" ] l|2008l ). Left: 
The curves are for our model with a = 3 and a range of values of Vcfi = 30, 50, 70 km . Right: The curves are for our model with 
''c.o = 50 km s~^ and a range of values of o = 1, 2, 3, oo. 



in Figure 9 of iRauch et alJ ()2008h . shown in black. The red 
dot-dashed curve shows the model assuming a duty cycle. 
The green dashed curve is for the reduced luminosity. 

So far we have assumed (following 
IPiikstra. Haiman. fc Spaand l2006l ) that all gravitational 
binding energy is converted into Lya radiation. In reality, 
hot accretion flows could result in substant ial amounts 
of b remsstrahlung and He"*" line emission (|Keres et al.l 
at the expense of I/Lya- Furthermore, the observed 
luminosities have not been corrected for slit losses which 
can be as large as a factor of a few. We therefore conclude 
that cooling radiation is indeed unlikely to contribute 
significantly to the majority of the faint Lyg emitters . 
This further supports the suggestion bv lRauch et"al] (|2008l ) 
that the emitters are predominantly powered by centrally 
concentrated star formation surrounded by extended Lya 
halos. By st a cking the spectra of a subset of the emitters, 
iRauch et al.l l|2008r ) showed that the emission from these 
halos appears to extend to radii even larger than those for 
the individually detected emission plotted in Fig. 2 by a 
factor of at least two. 



3.2 A simple model for the Lya luminosity 
function 

We will demonstrate now that a simple model where the Lya 
luminosity due to stars Lf^ya is proportional to the total 
mass {M oc v^) of the halos with virial velocity above a 
threshold fmin. 



-'Lya 



LI 




erg S ^ if Uc > Wmin 

otherwise. 



(15) 



fits the observed luminosity function remarkably well. Note 
that detailed numerical simulations of much brighter Lya 
break galaxies and Lya emitters appear to be consistent 
with this simple scaling of t he Lyg emission with t he prop- 
erties of the DM host halo (jNagamine et all 1200^ 1. In the 



right panel of Figure [31 we compare the observed lumi- 
nosity function of the faint Lya emitters with a luminos- 
ity function modelled in this way for a range of values of 
fd- The observed luminosity function is fit well for the fol- 
lowing parameter combinations (/d, Lq) '^min [ km s-i]) = 
(1,4 X 10'*\ 75), (0.5,8 X 10*S 60), (0.2, 1.8 x 10''^45). The 
upper and lower boundary of the grey shaded region can be 
fit with values of Lq and Vmin that differ from the quoted 
values by ~ 25% and ~ 5% respectively. 

Our assumed scaling I/Lya oc M is shallower than would 
be inferred from the star forma tion in most models of galaxy 
formation. The simulations of iPontzen et al.1 l|2008l ) would 
e.g. predict Lf^ya oc These models are, however, gen- 

erally tuned to produce a rather shal low faint end of th e 
luminosity function. As discussed by I Ranch et all ||2008D . 
their observations imply that the luminosity function of Lya 
emitters steepens cons iderably at very faint luminosities (cf. 
iLe DeUiou et al.ll2006l '). 

Intriguingly, the values for the velocity cut-off' and duty 
cycle in our model are very similar to those required to fit the 
velocity width distribution of the associated low-ionization 
metal absorption of DLAs and the size distribution inferred 
from the Lya emitters. We caution, however, that the sig- 
nificance of the apparent turn-over at ~ 1.25 x 10*° erg s~^ 
at the faint end of the observed luminosity function is un- 
certain. At faint ffux levels the luminosity function will be 
affected strongly by the sensitivity limit of the observations. 



4 THE MASSES AND VIRIAL VELOCITIES 
OF DLA/LL HOST GALAXIES 

If the population of faint Lya emitters detected by 
iRauch et all l|2008l ) can be identified with DLA/LL host 
galaxies, then this constitutes the first measurement of the 
space density and average size of DLA/LL h ost galaxies. 
The last section supported the suggestion by iRauch et al.l 
that the Lya emission is powered by star formation. In this 
case, with standard assumptions for the conversion of Lya 
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emission to star formation rate, the Lya luminosities corre- 
spond to star formati on rates of 0.07 - 1.5 M^yr"^ , simi - 
lar to that inferred bv lWolfe. Prochaska. fc Gawised (|2003l ) 
from the CII* A1335.7 absorption in DLAs. No continuum 
is detected, so there is no information about stellar or total 
masses of the object. Our modelling is thus currently the 
only handle we have on the masses (and virial velocities) of 
what should be a statistically representative sample of DLA 
host galaxies. 

In Figure [l] we show the contribution of DM halos of 
different masses and virial velocities to the incidence rate 
of DLAs in our model, for the range of parameters used 
to model the suppression of the cross-section for damped 
Lya absorption in low-mass DM halos. We also show the 
results fr om two recently published numerical simul ations 
of DLAs (|Razoumov et al.ll2007l : iPontzen et al] 120081 '). The 
differential line density for DLAs is calculated similarly to 
A*', except that we consider intervals of dX and dlogj^Q M, 



d^AA 



dX d login M 



= — InlO M nM{M,X) a{M,X). 
Ho 



(16) 



Note that larger values of a result in a sharper turnover 
at low masses, at the expense of increasing the abundance 
of DLAs with high masses (the area under the plot is nor- 
malised) . 

The majority of the DLAs in our model have virial ve- 
locities in the range 50 to 200 km s~^, corresponding to 
total masses of 10^" to W^^Mq. As discussed in the previ- 
ous sections, the turn-over at small virial velocities is most 
constrained by the velocity width distribution of the asso- 
ciated low-ionization metal absorption and is most likely 
attributable to feedback effects due to star formation. The 
decline at large virial velocities and masses is due to the 
decline of the space density of DM halos. 

The incidence rat e in the numerical simulations of 
iRazoumov et al] (|2007l ) shows a similar peak, albeit shifted 
to somewhat smaller masses/virial velocities than our model 
requires to fit the kinem atical data of th e DLAs. This is 
perhaps not surprising — IRazoumov et all (|2007l ) find that 
the velocity widths in their simulations fall somewhat short 
of those observed. Their simulation also takes into account 
DLAs that are not contained within any halo i.e. intergalac- 
tic DLAs. 

The numerical simulations of iPont zen et all (|2008l ) 
show a sharper peak centered on virial velocities of 30-80 
km s~^. Such a sharp peak at rather low virial velocities ap- 
pears, however, at odds with the observed velocity widths of 
the associ ated low-ionizat ion metal absorption. In the sim- 
ulations of IPontzen et al.l . the decline of the contribution to 
the incidence rate with increasing mass is much faster than 
the decline of the space density of massive halos. This fast 
decline is due to a flattening of the absorption cross sec- 
tion with increasing mass in massive halos ([Pontzen et al.l . 
Figure 4). 



5 SUMMARY AND CONCLUSIONS 

We have considered an updated version of the Haehnelt et 
al. model for the kinematics of DLAs, in light of the dis- 
covery of a new population of extended low surface bright- 
ness Lya emitters with a total inferred incidence rate similar 



to that of DLAs. The main differences with the modelling 
of iHaehnelt. Steinmetz. fc RauchI (HqoO) are the use of the 
Sheth-Tormen modification to the Press-Schechter formal- 
ism, an update of cosmological parameters, and the use of 
an exponential suppression of the cross section for damped 
absorption for low virial velocities instead of an sharp cut- 
off. 

Our main results are the following. 



• The observed velocity width distribution of the associ- 
ated metal absorption can be fit with a model where the 
cross section for damped absorption scales with the virial 
velocity of the halo as a oc v'^'^ , the absorption cross section 
is suppressed in halos with Vc ^ 70 km s~^, and the con- 
ditional velocity width is given by that in the simulations 
of Haehnelt e t al. (or the very si milar distribution of the 
simulations of IPontzen et al.ll2008l ). 

• The same model can fit the cumulative incidence rate 
for absorption inferred from the size distribution of the Lya 
emitters if the Lya emission has a duty cycle of /d = 0.2 — 
0.4, and the emission extends over an area which is larger 
than the cross section for damped Lya absorption by a factor 



/d 



2.5- 



• The maximum expected Lya cooling luminosity due to 
collapsing gas in DM halos falls short of the observed Lya 
luminosities by a factor of three to five, even for optimistic 
assumptions regarding the expected emission due to Lya 
cooling. Furthermore, the expected dependence of the Lya 
cooling luminosity on the virial velocity of DM halos thereby 
maps into a luminosity function with a slope shallower than 
observed. Lya cooling radiation should thus not contribute 
significantly to the Lya emission for the majority of the 
objects, especially at the faint end of the Lya luminosity 
function. 

• The Lya luminosity function is well fit by a simple 
model where the Lya luminosity scales linearly with the 
mass of the DM halo and the emission is suppressed for low 
mass DM halos. The Lya luminosity function can be fit for 
a wide range of duty cycles including the duty cycle required 
to simultaneously explain the kinematic properties of DLAs 
and the cumulative incidence rate inferred from the observed 
size distribution of the Lya emitters. 

• Our model predicts that the bulk of the contribution to 
the incidence rate of DLAs comes from absorption systems 
hosted in DM halos with virial velocities in the range from 
50-200 km s"^ and masses in the range 10^° to IO^^Mq. The 
cut-off for damped absorption occurs at somewhat higher 
virial velocities than suggested by numerical simulations, 
which attempt to simulate the gas distribution and kine- 
matics of DLAs self-consistently. These simulations, how- 
ever, fall short of reproducing the observed velocity widths 
of the associate low-ionization metal absorption of DLAs. If 
the suppression of the cross section for damped Lya absorp- 
tion in halos with virial velocities up to 70 km s~^ is indeed 
real, then feedback due to star formation at high redshift 
has to be more efficient in removing gas — even from rather 
deep potential wells — than is assumed in most models of 
galaxies formation and numerical simulations. Alternatively 
the simulation (and our model) may underestimate the ef- 
fect of stellar feedback on the velocity width of absorbers 
hosted by DM halos with small virial velocities. This may 
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not be implausible given the fact that the simulations still 
fail to reproduce realistic galactic winds. 

Our modelling here adds to the mounting evidence that 
DLAs and LLs predominantly probe the outer regions of 
dwarfislQ galaxies. This picture had been suggested for a 
while by the kinematic and metallicity data. With the dis- 
covery of a faint population of Lya emitters, most plausibly 
identified with the population of DLA/LL host galaxies, we 
finally have a handle on their space density, sizes and (in 
conjunction with models like the one presented here) masses 
and virial velocities. In the picture that emerges, DLAs are 
hosted by the galaxies that, in the context of the now well 
established ACDM paradigm for structure formation, are 
expected to become the building blocks of typical galaxies 
like our Milky Way. 
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* "In English the only correct plural of dwarf is dwarfs, and the 
adjective is dwarfish . . . dwarves and dwarvish are used . . . only 
when speaking of the ancient pe ople to whom T horin Oakenshield 
and his companions belonged. " iTolkienI lll966h . 



